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PREFACE 


The  Defense  Advanced  Research  Projects  Agency  (ARPA)  has  for  sev¬ 
eral  years  (through  1977)  supported  the  design  and  technical  evaluation 
of  a  passive  communications  satellite  (PACSAT)  by  the  Stanford  Research 
Institute  (now  SRI  International).  In  1977,  the  Defense  Communications 
Agency  (DCA)  asked  The  Rand  Corporation  to  examine  the  utility  of  PACSAT, 
as  well  as  other  passive  satellite  designs  to  support  a  communications 
link  between  the  National  Command  Authority  and  the  airborne  strategic 
bomber  force.  The  present  study  examines  the  use  of  PACSAT  in  a  com¬ 
munications  link  to  the  MX  mobile— missile  sites.  That  application 
differs  from  those  of  the  earlier  Rand  study  in  that  the  location  of 
the  MX  missile  sites  is  known  to  the  transmitter,  thus  allowing  a  priori 
carrier  frequency  selection  and  relaxing  the  antenna  requirements  for  an 
ideal  PACSAT  configuration.  In  addition,  this  study  examines  the  per¬ 
formance  of  PACSAT  for  nonideal  configurations. 

The  study  reported  in  this  Note  examines  the  possibility  of  .using 
signal  processing  techniques  for  correcting  the  PACSAT  return  signal  for 
near-field  effects  and  signal  acquisition.  The  results  should  be  of 
interest  to  analysts  associated  with  strategic  communications,  planning, 
and  operations. 

James  C.  Springett  and  Irving  S.  Reed,  the  author  of  the  appendix, 
are  consultants  to  The  Rand  Corporation. 
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I .  INTRODUCTION  AND  SUMMARY 


The  concept  of  a  linear  discrete  element  reflecting  array  in  earth 
orbit  as  a  method  for  providing  reliable  communications  between  highly 
selective  geographical  areas  has  been  proposed  by  Yater  [1]  and  subse¬ 
quently  developed  by  SRI  [2].  Such  a  passive  communication  satellite 
(PACSAT) ,  for  example,  may  serve  to  provide  command  communications  for 
portions  of  the  strategic  forces.  For  a  sufficiently  large  (long) 
PACSAT,  individual  sites  may  be  addressed  because  the  return  signal 
beamwidth  is  very  narrow.  Addressing  is  accomplished  by  carrier  fre¬ 
quency;  the  nominal  requirements  for  such  services  will  be  covered  in 
a  companion  Rand  Note  [3] . 

This  Note  is  primarily  concerned  with  two  problems  that  affect 
efficient  PACSAT  system  operation:  (1)  the  need  for  the  receiving 
terminal  to  search  for  the  carrier  frequency  being  employed,  and 
(2)  a  method  by  which  the  full  potential  of  the  array  may  be  realized 
when  the  receiving  terminal  operates  in  the  near-field  of  the  reflected 
electromagnetic  radiation.  The  first  problem  results  from  an  assumed 
inability  of  the  receiving  terminal  to  know  exactly  the  location  (in 
terms  of  incident  angle  to  the  array)  of  the  transmitter  (plus  some 
additional  uncertainty  as  to  the  possible  tilt  angle — relative  to  the 
ideal  vertical— of  the  array).  As  a  result,  the  receiver  must  search 
(over  a  limited  range)  for  the  frequency  that  the  transmitter  uses 
to  illuminate  the  receiver.  The  second  problem  results  from  the 
length  of  the  array  and  the  use  of  a  low  orbit  altitude.  Without 
some  form  of  compensation  processing,  the  individual  reflected  waves 
from  each  element  of  the  array  arrive  at  the  receiving  antenna  with 
sufficient  cumulative  phase  error  that  the  potential  gain  of  the  array 
is  significantly  degraded. 

References  to  the  main  body  of  the  text  may  be  found  on  p.  37. 

**The  transmitter  is  assumed  to  know  exactly  the  location  of  the 
receiver  and  all  other  essential  parameter  values. 
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A  solution  to  the  frequency  search  problem  is  .to  use  a  discrete 
set  of  frequencies  spaced  by  an  amount  sufficient  to  steer  the  reflected 
beam  by  one  beamwidth  between  any  two  adjacent  frequencies.  Thus,  the 
transmitter  selects  the  particular  frequency  that  places  the  receiver 
within  the  return  beam,  and  the  receiver  searches  over  a  subset  of  the 
frequencies  to  locate  that  being  employed. 

To  minimize  the  near— field  loss,  it  is  proposed  that  the  carrier  be 

modulated  by  a  pseudo-noise  (PN)  code  which  permits  subdivision  of  the 

array  return  into  spatial  segments  by  means  of  correlation  techniques 

at  the  receiver.  Since  each  spatial  segment  acts  as  a  short  array  for 

which  the  near-field  cumulative  phase  error  is  small,  the  mean  phase 

differences  between  the  segments  may  be  corrected  (removed) ,  and  the 

segment  components  combined  to  effect  the  return  from  the  total  array 

A 

with  an  acceptably  low  loss. 

The  following  sections  develop  the  various  concepts.  Section  II 
introduces  an  electrical  transfer  function  model  for  the  array  which 
is  useful  to  the  ensuing  developments.  Section  III  discusses  the  fre¬ 
quency  acquisition  process,  and  Section  TV  addresses  the  near-field 
correction  processor.  Finally,  Section  V  concludes  with  an  outline 
for  further  work. 

The  appendix  discusses  an  alternate  approach  for  a  frequency  de¬ 
tector. 

''Compensation  of  this  form  was  originally  suggested  by  Dr.  Irving 
S.  Reed,  Rand  Consultant. 
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II.  PACSAT  ARRAY  TRANSFER  FUNCTION  MODEL 


Figure  1  shows  the  basic  geometry  of  the  transmitter  (T)  and  re¬ 
ceiver  (R)  relative  to  the  bottom  (nearest)  element  of  the  array  which 
is  located  at  an  altitude  H  above  the  perpendicular  to  the  baseline 
TR.  The  distance  from  T  to  element  1  is  designated  P^,  and  the  dis¬ 
tance  from  T  to  the  nth  element  is  p  .  Similar  designations  are  made 
for  the  R  to  element  distances,  as  shown.  Angles  to  the  array  are  de¬ 
fined  as  indicated  on  Fig.  1.  In  total,  there  are  N  elements  in  the 
array,  and  each  element  is  spaced  by  a  distance  s. 

The  distances  are  related  to  the  altitude  and  the  angles  in  the 
following  manner. 


=  V  H2  +  H2  tan2  a 

(1) 

"\/(H  +  ns)2  +  H2  tan2  a 

(2) 

Combining  (1)  and  (2)  yields: 


P 


In 


(3) 


By  expanding  the  radical  in  a  series,  and  neglecting  the  third  order 
and  higher  terms  which  have  no  practical  effect  for  the  PACSAT,  (3) 
becomes : 


Similarly, 


In 


Pi  + 


Hns 


(ns) 

2Pi 


2n 


Hns  (ns) 2 

P2  2p2 


(4) 


1 


(5) 


^1n  sin 

(H  +  ns)  tan  Qn 


Fig.  1  —  PACSAT  array  basic  geometry 


P2  sin  (3  =  P2n  sin  /3n 
=  H  tan  /3  =  ( H  +  ns }  tan  /3  n 
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Now,  defining  p  =  p1  +  p2  and  pn  =  pln  +  P2n>  the  differential  dis 
tance  between  the  first  and  nth  elements  is: 


The  differential  propagation  time  delay  is  obtained  by  dividing 
(6)  by  the  propagation  velocity  of  light,  c,  viz 


By  design,  the  time  delay  4-  is  to  be  such  that,  when  the  array 
is  illuminated  by  a  frequency,  fQ, 

f  (t  4-  x  )  =  1  .  (2) 

o.l  2 

Noting  from  Fig.  1  that  H/ =  cos  a  and  H/P2  =  cos  3, 

sf 

f  (x.  +  x „)  =  — —  (cos  a  +  cos  3)  =  1  .  (9) 

o  1  2  c 


When  a  =  6  =  0  (end-fire  condition),  2sfQ/c  -  1,  or  f  -  c/2s.  Desig¬ 
nating  this  particular  frequency  by  the  symbol  f f ,  i.e.. 


(10) 


and  substituting  this  back  into  (9),  the  fundamental  frequency  steering 
relationship  for  the  array  is  obtained: 


f 

o 


cos  a  +  cos  6 


(ID 
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The  selection  of  is  the  option  of  the  system  designer;  however, 
once  it  is  fixed,  then  the  inter-element  spacing  of  the  array  elements, 
s,  is  established  by  (10). 

Using  the  relationships  (8),  (10),  and  (11),  (7)  may  be  rewritten 

as : 

a  -  _EL  .  n2e  cos(a)  sinz(a/2)  +  cos(fi)  sin2(3/2) 

Tn  f  2f0f  H  cos(a)  +  cos(3) 

o  Z  o  L 

The  differential  phase  delay  of  the  carrier  at  frequency  f  =  / 2ir 

is  given  by  co  Ax  =  A <t>  ,  thus: 

J  o  n  n 

.  ,  „  ,  7Tn2c  Tcos(a)  sin2(a/2)  +  cos(3)  sin2(0/2) 

~  2™+-q*  If- - cos(a)  '+  c'osTO 

The  second  term  in  (13)  represents  the  quadratic  curvature  of  the  wave 
front  that  builds  along  the  array.  Differentially  between  successive 
elements  it  is  very  small;  however,  when  n  is  large,  its  value  rela¬ 
tive  to  the  first  element  cannot  be  ignored  for  certain  altitudes  H 
(especially  when  H  is  one  earth  radii  or  less) .  It  thus  manifests 
itself  as  the  near-field  problem.  For  the  sake  of  the  ensuing  dis¬ 
cussions,  it  will  be  given  the  symbol  <j>  ,  i.e., 

,  A  TTn2c  cos(a)  sin2(a/2)  +  cos(3)  sin^(3/2) 

^n  f^H  cos (a)  +  cos (3) 

Equations  (7)  and  (13) ,  together  with  the  connecting  relation¬ 
ships,  suggest  a  useful  signal  model  for  the  array,  shown  in  Fig.  2. 

* 

The  model  consists  of  a  pair  of  tapped  delay  lines  interconnected  by 
phase  shifters.  The  delay  sections,  and  x are  as  defined  in  (7), 
and  respectively  represent  the  inter-element  arrival  delay  time  for  a 
plane  wave  incident  from  the  angular  direction  a,  and  the  inter-element 
scattering  delay  time  for  a  plane  wave  received  (in  the  far-field)  at 

•k 

Finite  dimension  and  mutual  coupling  effects  between  the  elements 
have  been  ignored. 
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Fig.  2  —  Array  signal  transfer  function  model 


an  angular  direction  0.  At  the  carrier  frequency  fQ,  +  x2  effec 
tively  delays  the  wave  by  one  full  cycle  or  shifts  it  in  phase  by  2tt. 

The  model  uses  delay  sections  rather  than  phase  shifters  because,  in  gen¬ 
eral,  the  carrier  is  modulated  and  phase  shifters  alone  (valid  only 
at  f  )  are  insufficient.  On  the  other  hand,  the  quadratic  delay  com¬ 
ponents  (T  of  (7))  are  extremely  small  and  do  not  represent  any  real 
n  _  . 

effect  on  the  modulation.  Therefore,  it  is  appropriate  to  model  the 

x  as  phase  shifts  qnly  of  the  carrier  frequency,  each  phase  shift, 
n 

4>  ,  being  given  by  (14)  . 

If  a  signal  S  (t)  =  m(t)  cos^t)  is  input  to  the  array  model  of 
Fig.  2,  the  received  signal  is  given  by: 


N 

S  (t)  =  V  m[t  -  (n  -  1)t]  cos  [oj  t  -  c|>  ] 

R  M\«J 

n=l 


(15) 


where  T  =  T]_  +  x2*  t^ie  modulation  m(t)  varies  slowly  so  that  its 

delay,  for  n  ^  N,  can  be  ignored,  then  (15)  simplifies  to: 


N 

SR(t)  =  m(t)  2  cos  [WQt  -  • 

n=l 


(16) 


Finally,  for  an  ideal  array  operating  under  far-field  conditions,  the 

(j)  are  negligible  (<j>  =  0)  and  the  received  signal  becomes: 

n  n 


N 

S  (t)  *=  m(t)  cos(&)  t)  =  Nm(t)cos(0J  t) 

n=l 


(17) 


The  right-hand  expression  of  (17)  represents  perfect  in-phase  coherent 
combining  of  the  reflections  from  each  array  element. 

The  model  of  Fig.  2  with  all  ^  =  0  is  useful  for  deriving  cer¬ 
tain  properties  of  the  array..  The  impulse  response  of  the  array  for 
this  case  is  easily  written  by  inspection  as: 
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N 

h(t)  =  <5[t  -  (n  -  l)l],  (18) 

n=l 


where  5(t)  is  the  usual  Dirac  delta  function.  Taking  the  Fourier 
transform  of  (18)  produces  the  transfer  function  for  the  array,  viz., 

N 

H(oj)  =  f{  6  (t)  }  =  exp  [-  j  (n  -  1)tcoJ  .  (19) 

n=l 


Since  the  series  in  (19)  is  a  geometrical  progression,  it  is  easily 
summed  to  yield: 


H(oj) 


exp[-  jNTOj]  -  1 
exp  [  -  j  toj  ]  -  1 


The  magnitude-squared  of  H(w)  is  calculated  to  be  (w  -  2Trf)  : 


(20) 


|H(f)|2 


sin2(7TNTf  ) 
sin2  (inf  ) 


(21) 


Clearly,  j H (oo)  j 2  is  a  periodic  function  of  f.  Since  the  context  calls 
for  operation  of  the  array  at  a  frequency  in  the  vicinity  of  1/x,  (21) 
can  be  rewritten  in  the  form: 


sin2  [ttNt (f  -  fQ)] 
'  “  sin2  [ttt  (f  -  F7] 


(22) 


where  f  =  1/t,  and  has  already  been  given  in  terms  of  the  end  fire 
frequency  by  (11) . 

As  one  example  of  the  utility  of  (22),  the  -3  dB  bandwidth  of  the 
array  is  calculated.  When  f  =  fQ,  | H ( f ) ]  =  N.  Thus,  it  is  necessary 
to  determine  the  values  of  f  for  which  i H ( f ) |  =  N//2  ,  i.e. ,  there 
exists  f,  <  f  and  fn  >  f  for  which 
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sin[7TNT(fo  -  fx)]  sin[uNT  (f^  -  fQ)]  N 

sin[TTT (f ^  -  f^)]  sinfTTTCf^  -  fQ)]  ^2 


(23) 


Letting,  say,  x  =  Trx(f  —  f^)  ,  it  is  sufficient  to  solve  for  x  in  the 
expression 


sin(Nx-)  _  N  -  (24) 

sin(x)  /2 

Because  the  bandwidth  is  narrow  compared  to  the  center  frequency  f q , 
x  will  be  small  (x  <<  it),  so  that  sin(x)  —  x,  and  (24)  can  be  rewritten 

as : 


sin(Nx)  _  _J_ 
(Nx)  /2 


(25) 


Reference  to  a  table  of  sin(z)/z  quickly  reveals  that  Nx  -  1.39.  Thus, 
Nttt  (fQ  -  fx)  =  1.39,  or 


f 

o 


1.39  =  1.39 
Ntit  Ntt  o 


(26) 


Now  f  -  f  is  one-half  of  the  -3  dB  bandwidth;  therefore,  the  total 
’  o  1 

-3  dB  bandwidth  is  twice  that  given  by  (26),  or: 


BW (-3  dB)  = 


2.78 

Ntt 


f 

o 


0.885 
N  o 


(27) 


Substituting  (11)  for  f  in  (27)  yields: 


BW(-3  dB)  = 


1-77  fl 

N(cos  a  +  cos  0) 


(28) 


This  result  shows  that  the  bandwidth  is  inversely  proportional  to  the 
length  of  the  array,  and  increases  as  the  angles  to  the  array  increase. 

Further  applications  of  the  model  and  relationships  developed  in 
this  section  will  be  made  in  Sections  III  and  IV. 
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•  III.  SIGNAL  FREQUENCY  ACQUISITION 

The  problem  of  frequency  acquisition  may  be  readily  appreciated 
by  resort  to  a  typical  use  scenario. 

An  airborne  command  post  (ABNCP)  is  to  communicate  sequentially 
with  M  MX  receivers  via  the  PACSAT  in  equatorial  orbit.  Depending 
upon  the  location  of  the  ABNCP,  the  MX  field,  and  the  PACSAT,  there 
exists  a  range  of  angles  a  and  3  which  encompasses  all  communication 
scenarios . 

At  any  instant  of  time,  the  ABNCP  knows  the  a  and  3  necessary  to 
place  a  given  MX  receiver  at  the  center  of  the  PACSAT  return  illumina¬ 
tion  annulus.  This  is  by  virtue  of  a  priori  MX  movements,  computer 
projections  of  the  PACSAT  location,  and  accurate  ABNCP  coordinate  deter 
mination.  The  MX  receiver  is  also  assumed  to  know  the  requisite  6, 
however,  a  is  generally  unknown  to  the  receiver  because  there  is  no 
realistic  way  to  communicate  this  information  from  the  ABNCP  to  the  MX 
site.  Thus,  the  receiver  must  expect  the  a  to  fall  on  some  range  of 
values  dependent  upon  the  immediate  3  and  the  limits  imposed  upon  ot 
from  flight  range  constraints  on  the  ABNCP  relative  to  some  geo¬ 
graphic  center. 

For  a  given  3,  there  is  a  corresponding  ct,  which,  from  the  re¬ 
ceiver’s  viewpoint,  falls  on  a  range: 

^  a  ^  a2  . 

Furthermore,  considering  the  total  MX  field  and  all  useful  PACSAT 
orbit  positions  (or  the  range  of  orbit  positions)  ,  3  could  be  any 
angle  on  the  range: 


3  .  ^3^3 

min  max 

The  a  uncertainty  translates  into  a  transmitted  frequency  uncer¬ 
tainty  at  the  receiver.  This  frequency  uncertainty  is  obtained  using 
(11)  as:  (a2  >  c^) 
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_ _ _ -  _  -  -  2f  (29) 

cos  +  cos  3  cos  +  cos  B  “  SL 

where  is  the  transmission  frequency  used  for  ’the  PACSAT  array  end- 
fire  condition  of  a  =  3  =  0.  Thus,  the  receiver  must  acquire  over 
this  range  of  frequency  uncertainty.  To  be  realistic,  the  frequency 
acquisition  should  not  significantly  extend  the  time  necessary  to 
transmit  the  basic  message  from  the  ABNCP  to  the  MX  receiver. 

The  proposed  solution  to  the  frequency  uncertainty  and  acquisi¬ 
tion  problem  is  to  use  a  fixed  set  of  selectable  frequencies  (at  the 
transmitter)  structured  in  the  following  manner.  Each  frequency  cor¬ 
responds  to  a  specific  ot  and  3^  pair.  The  frequencies  are  spaced  so 
that  for  any  exact  a  and  3,  there  is  a  frequency  f  (corresponding  to 

a  and  3  )  for  which  the  receiver  is  guaranteed  to  fall  within  the 
n  *  n 

-3  dB  beam  angle  of  the  return  radiation  from  the  array.  Stepping 

(by  the  transmitter)  from  one  frequency  to  an  adjacent  frequency  in 

* 

the  set  causes  the  return  beam  to  step  in  angle  by  one  -3  dB  beam- 

width.  Thus,  when  the  ABNCP  determines  the  a  and  3  to  communicate  with 

a  particular  MX  receiver,  it  selects  the  closest  and  3^  pair,  and 

uses  the  corresponding  transmission  frequency,  f  . 

The  receiver  is  assumed  to  know  3  (by  virtue  of  computer  prediction, 

or  direct  measurement  using  self-illumination  by  a  signal  transmitted 

to  the  array  from  the  receiving  site).  However,  only  the  nominal  a 

(ABNCP  geographic  center)  is  known  to  the  receiver,  and  there  is  an 

uncertainty  range,  0^  to  a  .  Since  the  receiver  knows  the  frequency 

set  and  the  3  used  by  the  transmitter,  it  need  then  only  search  over 
n 

the  candidate  f  corresponding  to  the  a  values  that  fall  within  the 
n  n 

uncertainty  range. 

To  accomplish  the  frequency  search,  it  is  postulated  that  the 
receiver  will  use  an  N^  channel  frequency-selective  energy  detec¬ 
tor.  Each  channel  is  set  to  an  expected  frequency,  and  has  band¬ 
width  sufficient  to  encompass  all  frequency  uncertainty.  (The  bandwidth 

Or  perhaps,  in  order  to  allow  for  some  margin,  the  -2  dB  beam- 
width  is  more  practical. 
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will  be  found  to  be  considerably  less  than  the  inter frequency  spacing.) 

The  integration  period  and  detection  threshold  are  selected  to  give  a 

high  detection  probability  (say,  PD  =  0.999)  when  a  signal  is  present 

in  one  of  the  channels,  and  a  very  low  false  alarm  probability  (say, 

-14 

P  =10  )  when  no  signal  is  being  received. 

FA 

This,  then,  is  the  broad  concept  for  frequency  acquisition.  It 
now  suffices  to  explore  some  of  the  details,  and  to  conclude  with  a 
numerical  example  to  illustrate  the  probable  size  of  the  system. 

To  begin,  the  sensitivity  of  the  return  beamwidth  as  a  function 
of  a  and  3  is  derived.  The  array  illumination  frequency,  fQ,  is  given 
by  (11)  for  a  specified  a  and  3.  It  should  be  recalled  that  when  (11) 
is  satisfied,  the  inter-element  propagation  times  for  a  plane  wave 
arriving  from  angular  direction  a,  and  culminating  in  a  far-field 
plane  wave  being  received  in  the  angular  direction  3,  is  exactly  one 
cycle,  or  a  2ir  phase  delay  of  f  .  Now,  if  the  receiver  is  not  located 
exactly  at  the  angle  3,  but  rather  at  an  angle  3  +  A3 ,  the  phase  shift 
between  the  individual  returns  from  each  element  will  no  longer  be 
exactly  2TT,  and  a  reduced  resultant  combined  amplitude  will  be  ob¬ 
tained.  In  fact,  for  certain  A3,  nulls  (no  resultant  signal)  are 
manifest. 

The  beamwidth  may  be  derived  by  starting  with  the  relationship 
(obtained  from  the  definition  of  T ^  in  (7),  plus  (10)  and  cos  3  =  H/p^)  : 

cos  3  =  2f£T2  .  (30) 

By  deduction, 

cos(3  +  A3)  =  2fji(T2  “  =  cos  ^  “  2f£Ai2  ,  (31) 

and  solving  for  A3  yields: 

A3  =  cos  ^  [cos  3  -  2f^Ax2]  -  3 
If  the  beamwidth  is  defined  by  63  =  2A3 ,  then: 


(32) 
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53  =  2  cos  ^  [cos  3  -  -  23  .  (33) 

It  now  remains  to  derive  At^.  From  the  array  model  of  Fig.  2,  with 

the  (p  '  s  =  0  and  T0  replaced  by  T  +  At  ,  the  received  signal  is 
n  2  2  2 

written  as: 

N 

S  ft)  -  Y  cos [oj  t  +  noj  (t.  +  T0  +  At)]  .  (34) 

R  Z— t  o  o  1  2 

n=l 


Recalling  that,  by  design,  C0o(t^  +  T?)  ~  2tt ,  (34)  reduces  to: 


N 

S  (t)  =  co s [oj  t  +  nu  At.] 

R  Lu  o  o  2 

n=l 

N 

=  cos  Oi  t  co s ( noo  At„) 

o  l—i  o  2 

n=l 

N 

-  sin  (o  t  )  sin(nto  At  ) 
n=l 


=  A  cos  oj  t 
c  o 


A  sin  oo  t  . 
s  o 


(35) 


Now,  when  Ato  =  0,  A  =  N,  and  A 
*  2  c  s 


0.  Thus,  the  condition  for  which 


V 


c 


+  A2 
s 


_N 

/2 


(36) 


must  be  established.  Summing  the  series  representing  k^  and  Ag , 
may  be  readily  shown  that : 


it 
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+  A2 
s 


/Noi  At, 


o  2 

sin 

2  } 

NO)  At 
o  2 


(37) 


Equating  (36)  and  (37),  and  solving  for  At^  gives 


At 


2 


2.78 

Noj 

o 


(38) 


Taking  this  result,  together  with  (11)  and  (33),  produces  the  final 
result  for  the  beamwidth,  viz.. 


63  -  2  cos 


-1 


cos 


1.39 

UN 


(cos  a  +  cos  8) 


-  28 


(39) 


A  study  of  (39)  discloses  that  58  is  maximum  for  a  =  3  =  0  (end-fire) 
and  decreases  for  increasing  a  and  8.  Since  for  the  PACSAT  scenario 
the  end-fire  condition  does  not  usually  occur,  a  more  utilitarian 
approximation  to  (39)  may  be  derived.  When  8  >>  A3,  (31)  may  be 
written  approximately  as: 


or 


cos  8  -  A8  sin  3  =  nos  8  -  2f ^At ^  , 


(40) 


A8 


2££At2 
sin  8 


(41) 


Using  (11)  and  (38)  together  with  the  definition  for  63,  the  approxi¬ 
mate  and  simplified  result  is  obtained  as: 

?-q  ~  2 . 78  [cos  a  +  cos  3 
“  UN  L  sin  3 


3  >>  63 


(42) 


(thus  verifying  Eq.  (10)  of  [2]). 
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Recalling  from  the  earlier  discussion  that  steering  the  beam 
by  means  of  a  discrete  set  of  transmitter  frequencies  is  the  goal, 
it  is  now  of  interest  to  determine  the  spacing  and  number  of  frequen¬ 
cies.  A  step  from  one  frequency  to  the  next  will  be  assumed  to  step 
the  beam  in  angle  by  From  (11) , 


f  +  6f 

o 


_ 2JjL _ 

cos  a  +  cos(3  +  53) 


(43) 


2f 


£ 


cos  a  +  cos  3  -  53  sin  3 


(44) 


Using  (42)  to  substitute  for  63  sin  3,  and  simplifying,  the  approxi¬ 
mate  result  becomes: 


6f 


2.78 

ttN 


2f , 


cos  a  +  cos 


(45) 


Again,  recalling  that  the  transmitter  frequency  adjustment  range, 
needed  to  place  the  receiver  within  the  -3  dB  beamwidth,  is  Af  given 
by  (29),  the  number  of  discrete  frequencies  required  becomes: 


Nf  = 


Af 

5f 


7TN 

2.78 


(cos  a  +  cos  3) 


1 


cos  +  cos 


cos  a,  +  cos  3 


(46) 


An  example  at  this  juncture  will  serve  to  illustrate  the  results 
of  Eqs.  (45)  and  (46) .  A  satellite  1500  meters  in  length,  designed 
to  operate  at  an  end-fire  frequency  of  8  GHz,  will  have  N  =  80,001 
elements.  For  a  one-earth-radius  equatorial  orbit  altitude,  typical 
a  and  3  values  are  a  =  22°  and  3  =  35°  (using  the  ABNCP  and  MX  scenario). 
Assuming  =  20°  and  -  24°,  the  frequency  set  frequency  spacing  is 
found  to  be  6f  =  101.3  kHz,  and  the  number  of  frequencies  is  calculated 
as  N^.  =  1355. 

Attention  is  now  turned  to  the  means  by  which  the  MX  receiver  may 
be  brought  to  a  message  receiving  state. 
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The  frequency  acquisition  process  by  the  MX  receiver  consists  of 
two  steps:  (a)  detection  of  which  frequency  the  transmitter  is  using 
and  (b)  phase-locking  on  the  detected  frequency. 

Figure  3  shows  a  configuration  for  the  frequency  detector.  A 
channel  is  provided  for  each  of  the  candidate  frequencies,  and  each 
channel  is  functionally  identical. 

The  accuracy  of  any  received  frequency,  relative  to  its  theoreti¬ 
cal  static-condition  (time-instant)  value,  depends  on  the  transmitter 
and  receiver  frequency  stabilities  and  baseline  accuracies,  as  well  as 
the  degree  of  doppler  compensation  at  the  transmitter.  Generally, 
the  input  bandpass  filter  (BPF)  to  each  channel  of  the  frequency  de¬ 
tector  has  a  bandwidth,  W,  with  property 

W  «  Sf,  (47) 


as  W  need  only  be  sufficient  to  accommodate  the  frequency  uncertainty. 
Threshold  TH  is  selected  to  obtain  the  desired  false-alarm  probability 
(P  )  in  the  absence  of  a  signal  within  a  channel,  and  the  required  de¬ 
tection  probability  (P  )  when  a  signal  is  present.  For  the  absence 
of  signal,  all  y  outputs  take  on  value  -1,  and  when  a  signal  is  pre¬ 
sent  in  a  channel,  that  y  output  has  value  +1. 

The  integration  time,  T,  is  the  dependent  quantity,  and  is  re¬ 
lated  to  the  other  governing  parameters  by  the  relationship  [4] : 

T  =  Wd2/(C/N  )2  ,  (48) 

o 

where  d  =  Q  (P^)  -  Q  (P^)  ,  (49) 


and 


C/N  =  the  received  signal-to-noise 
spectral  density  power  ratio. 


Q  1(  )  is  the  inverse  normal  cumulative  distribution  function. 

For  a  practical  acquisition/detection  procedure,  an  overhead  of 
3T  seconds  should  be  allowed  for  detection.  T  seconds  are  needed 

"The  appendix  investigates  an  alternate  approach  that  is  simpler  to 
implement  but  requires  significantly  increased  decision  time  compared 
with  tha  proposed  configuration. 


Fig.  3  —  Functional  frequency  detector 
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for  detection,  T  seconds  for  detection  confirmation  (it  is  assumed 
that  two  successive  U  -  1  conditions  from  the  same  channel  are  neces¬ 
sary  to  proceed  to  the  next  step) ,  and  T  seconds  for  lack  of  signal 
onset  epoch  synchronization. 

Once  detection  has  been  accomplished,  a  separate  Costas— loop  type 
receiver  [5]  is  frequency  programmed  to  the  nominal  channel  frequency. 

This  loop  must  then  acquire  the  "carrier"  from  the  range  of  frequency 
uncertainty.  (The  word  carrier  is  placed  in  quotation  marks  because 
the  received  signal  is  assumed  to  be  of  a  suppressed  carrier  form, 
and  the  Costas-loop  regenerates  a  double-frequency  carrier  through 
the  inherent  squaring  process.)  Generally,  it  may  be  assumed  that 
the  frequency  uncertainty  will  be  greater  than  the  data  (or  symbol) 
rate  that  modulates  the  carrier.  As  a  result,  Costa$-loop  false  lock 
at  frequency  offsets  of  multiples  of  one— half  of  the  data  rate  may  be 
expected.  To  preclude  such,  it  is  expedient  for  the  transmitter  to 
remove  all  data  modulation  and  transmit  an  unmodulated  carrier  during 
both  the  frequency  detection  and  Costas-loop  acquisition  phases  of 
the  reception  process.  Thus,  unmodulated  carrier  is  transmitted  by 
the  ABNCP  for  a  prescribed  period  to  be  defined  subsequently. 

The  Costas  loop  may  be  acquired  using  one  of  two  basic  approaches  [6]: 
(a)  self-pull-in  and  (b)  frequency  sweeping.  The  former  is  the  sim¬ 
plest  to  implement  as  nothing  is  required  outside  of  the  loop  itself. 
However,  the  price  paid  is  generally  a  moderate  pull-in  time.  On  the 
other  hand,  the  sweep  approach  results  in  relatively  rapid  acquisition, 
but  at  increased  complexity  due  to  the  addition  of  sweep— related  cir¬ 
cuits.  The  maximum  frequency  uncertainty  is  equal  to  one-half  of  the 


frequency  detector  channel  bandwidth,  i.e.,  f^  =  W/2.  For  a  self¬ 
pull-in  Costas-loop  operating  at  a  high  loop  SNR  (say,  15  dB) ,  the 
necessary  conditions  are  (the  loop  damping  factor  is  taken  as  unity): 


(51) 
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where  K  is  the  open-loop  gain  and  2B^  is  the  loop  two-sided  noise 
bandwidth  in  Hz. 

For  the  frequency  sweeping  scheme,  the  sweep  rate  f  is  given 

by: 


f  £0.112  (2BT)2  Hz/sec.  (52) 

sw  L 

Equality  results  in  a  0.9  probability  of  acquisition. 

If  the  sweep  profile  is  taken  as  a  linear  frequency  sweep  from 
one  edge  of  the  band  W  to  the  other,  and  then  a  return  to  the  center 
of  the  band  (i.e.,  removal  of  the  sweep  for  the  nominal  frequency 
value),  the  acquisition  time  is  given  by: 

T  =  1.5  W/f  .  (53) 

acq  sw 

To  obtain  a  measure  of  the  relative  time  involved  for  both  ac¬ 
quisition  methods,  the  ratio  of  (51)  to  (53)  is  taken  for  the  maximum 
allowable  sweep  rate  [equality  in  (52)],  thus: 


T  (self-pull-in) 
acq _ _________ 

T  (swept  loop) 
acq 


1.44  W 
2Bt 


(54) 


Clearly,  whenever  W  is  significantly  larger  than  2B^,  the  swept  ac¬ 
quisition  approach  is  essential  to  acquisition  time  minimization. 

For  the  PACSAT  applications  under  consideration,  W  is  expected  to  be 
at  least  an  order  of  magnitude  larger  than  2B^.  Swept  frequency  ac¬ 
quisition  will  therefore  be  assumed. 

To  close  this  section,  consideration  is  given  to  the  total  time 
required  by  the  ABNCP  to  communicate  an  bit  message  to  a  selected 
MX  site.  The  total  time  is  that  necessary  for  frequency  detection, 
frequency  acquisition,  symbol  synchronization,  and  the  actual  message 
content.  Frequency  detection  and  acquisition  time  have  already  been 
discussed. 
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Symbol  synchronization  and  message  times  are  measured  in  bits. 

If  an  alternating  +1  and  -1  pattern  is  used  to  rapidly  establish 
symbol  synchronization  by  means  of  a  digital  transition  tracking 
loop  (DTTL)  synchronizer,  the  maximum  time  necessary  to  acquire  is 
about  120  symbol  periods.  With  a  symbol  rate  of  Rg  symbols  per  second, 
then  120/ Rs  is  the  absolute  time  necessary  to  acquire  symbol  synchroni¬ 
zation. 

Figure  4  shows  the  time  components  necessary  to  transmit  and 
receive  one  total  message.  In  order  to  virtually  guarantee  the 
Costas-loop  swept  frequency  acquisition,  the  maximum  allowable  sweep 
rate  [equality  in  Eq.  (52)]  has  been  backed  off  by  20  percent. 

At  this  point,  a  numerical  example  again  illuminates  the  reality 
of  the  total  communication  time.  The  following  assumed  parameter 
values  are  considered  typical  of  the  ABNCP/MX  scenario. 


N  =  400  bits 
M 

Rg  =  1000  symbols/sec 

C/N  =  40  dB  -  Hz 
o 

W  =4  KHz 


2B  =  250  Hz 

Li 


D 


=  0.999 
-14 


=  10 

FA 


d  =  10.75 


With  reference  to  Fig.  4,  the  total  time  is  calculated  to  be 
1.634  sec,  composed  of  the  following  components: 


Frequency  detection  time  =  0.014  sec 
Costas-loop  acquisition  time  =  1.1  sec 
Symbol  synchronization  time  =  0.12  sec 
Message  time  =  0.4  sec. 

It  should  be  noted  that  the  largest  component  is  the  Costas-loop  ac¬ 
quisition  time.  Practically,  this  can  only  be  made  smaller  by  taking 
steps  to  reduce  frequency  uncertainty. 
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Detection 

Acquisition 

Synchronization 
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"  *  ~ 

3Wd  2  /(C/N0)2 

16.7  W/(2Bl)2 

120/RS 

No  Symbol  Alternating 

Modulation  Symbols 


Fig.  4  —  Total  message  time  components 


As  a  final  observation,  it  is  seen  from  the  above  that  the 
smallest  component  of  the  total  communication  time  is  that  for  fre¬ 
quency  detection.  At  a  small  expansion  in  total  time — from  1.634  sec 
to  1.76  sec — the  frequency  detection  time  is  increased  by  a  factor 
of  10.  Through  this  modest  increase,  an  effective  hardware  savings 
is  obtained  in  the  number  of  physical  detection  channels  required  by 
the  frequency  detector  from  1355  channels  (see  earlier  example)  to 
49  channels.  Frequency  detection  is  then  effected  by  hopping  the  49 
channels  over  the  total  set  of  1355  frequencies. 


IV.  PROCESSING  FOR  NEAR-FIELD  PHASING  ERROR 


When  the  PACSAT  is  in  low  earth  orbit  (one  earth  radius  or  less) , 
and  the  array  is  sufficiently  long  (say,  1500  meters),  the  quadratic 
curvature  of  the  wavefront  that  builds  across  the  array  becomes  signif¬ 
icant  to  the  extent  that  the  receiver  is  not  in  the  far-field  of  the 
reflected  wave.  Thus,  unless  some  means  is  provided  to  correct  the  in¬ 
cremental  misphasing  of  the  return  waves  from  each  element  of  the  array, 
an  effective  loss  of  signal  power  results  from  the  quadratic  terms. 

A  second  source  of  phasing  error  results  from  a  mechanically  non¬ 
ideal  array.  Because  of  various  forces  acting  on  the  array,  the 
individual  elements  are  expected  to  be  displaced  (with  constraints) 
from  their  ideal  locations,  as  shown  in  Fig.  5.  Although  these  dis¬ 
placements  will  be  small,  they  could  represent  equivalent  phase  errors 
of  many  degrees.  Furthermore,  the  displacements  must  be  treated  as 
random  (although  they  will  be  correlated  over  short  segments  of  the 
array)  and  time  varying. 

Before  proceeding  with  the  development  of  a  solution  to  the 
phasing  problem,  an  example  of  the  magnitude  of  the  near-field  loss 
is  illustrative.  Equation  (14)  is  rewritten  in  the  form: 

=  £n2  (55) 

n 

with 

c  A  ttc  fcosCoQ  sin2(a/2)  +  cos(B)  sin2(3/2) 

£  f^H  |_  cos(a)  +  cos(3) 

The  coefficient  £  increases  with  increasing  a  and  3.  For  the  sake 

of  the  example,  a  and  3  are  both  taken  to  be  30°.  Additionally, 

H  =  1  earth  radius  (6378  km)  and  f.  =  8  GHz.  Calculation  of  £  using 

-9  ^ 

(56)  gives  £  =  1.237  x  10  radian.  The  effective  signal  loss  is  now 
calculated  using  the  array  center  element  as  reference  for  the  nominal 
a  and  3  angles,  yielding: 


(56) 
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Element  ideal  location 


Element  real  location 


Fig.  5  —  Nonideal  array  having  displaced  elements 
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(N-D/2 

1  +  2  ~y  ]  cos(en2) 

LoSS  =  - -  •  qi - - - —  ,  (57) 

N 

where  N  is  the  number  of  array  elements.  Taking  N  =  80001  (1500  meter 
array  length),  the  loss  obtained  from  (57)  is  -3.4  dB .  Clearly,  this 
loss  is  of  sufficient  magnitude  to  justify  the  use  of  a  corrective 
procedure. 

The  preceding  example  assumed  continuous  wave  (cw)  illumination  of 
the  array.  In  order  to  simply  develop  the  near-field  processing  con¬ 
cept,  assume  that  in  place  of  a  continuous  RF  wave  a  pulsed  RF  wave 
is  used,  with  a  pulse  width  equal  to  the  propagation  delay  time  T ^  + 
[given  by  (9)],  and  a  pulse  repetition  frequency  (PRF)  equal 
to  the  inverse  of  the  propagation  delay  for  the  entire  array,  i.e., 

PRF  =  [ (N-l)  ( T  +  t2)]_1  .  (58) 

With  such  pulsed  illumination,  the  received  wave  is  composed  of  a 
series  of  phase-disjoint  sequential  reflections  from  the  individual 
elements,  as  illustrated  in  Fig.  6.  It  is  seen  that  the  received 
signal  is  thus  composed  of  segmented  terms,  each  with  increasing  phase 
error  in  accord  with  Eq.  (14). 

Figure  7  shows  a  conceptual  near— field  processing  receiver  that 

corrects  for  the  components.  A  voltage-controlled  phase  shifter 
n 

is  used  to  sequentially  shift  each  receiv-ed  wave  segment  by  -(J)^  + 
so  that  at  the  output  a  continuous  phase  signal  (with  reference  phase 
(J)  )  is  obtained.  Thus,  all  of  the  deterministic  phase  error  is  effec¬ 
tively  removed.  (Note:  random  phase  errors  cannot  be  processed  away 
by  this  simple  approach  because  they  are  not  known  a  priori  at  the 
receiver. ) 

Although  the  foregoing  approach  to  solving  the  near-field  problem 
works  in  principle,  it  is  tar  from  practical.  Consider  the  following 
requirements  based  upon  the  previously  cited  system  parameters  (a  = 

8  =  30°,  N  =  80001,  S  =  1.875  cm).  The  pulse  width  is  calculated  to 


0  (t^  +  t2)  2<t1  +  t2)  3(t,  +  t2)  4  (Tj  +  r2 ) 


I 

r-o 

O' 

l 


Fig.  6  —  Pulsed  RF  wave  propagation  states  and  received  wave  structure 


-28- 


be  Ti  +  t2  =  °*1  ns’  and  PRF  =  115*5  KHz"  Also>  if  each  Pulse  is  to 
suppiy  the  same  received  energy  as  that  of  the  cw  signal  naving  power 

P  ,  then  the  pulse  power  is  given  by: 
cw 


P 


pulse 


(N  -  1)P 

cw 


(59) 


Typical  communication  link  calculations  show  that  P^  should  be  on 
'the  order  of  100  W.  Equation  (59)  then  requires  that  Ppulse  =  8  MW ! 
Clearly,  the  required  pulse  width  and  power,  together  with  the  im¬ 
plied  switching  and  processing  speeds,  render  the  pulsed  wave  concept 
unachievable. 

Before  proceedirig  with  a  description  of  a  practical  near— field 
correction  method,  the  requirement  for  correction  is  reevaluated. 

The  value  of  e  =  1.237  x  10~9  radian  is  very  small,  so  small,  in  fact, 
that  correction  of  phase  on  an  element— to— element  basis  is  unwarranted. 
Therefore,  it  is  reasonable  to  view  the  total  array  as  a  spatial  series 
of  subarrays,  each  having  N'  elements.  The  combined  returns  from  each 
subarray  are  then  accepted  with  some  loss,  and  near— field  phase  cor¬ 
rection  is  made  only  to  the  mean  phase  of  the  return  from  each  individual 
subarray.  In  order  to  assess  how  large  a  subarray  may  be,  consider 
the  following  table  of  near-field  dephasing  losses. 


N'  Dephasing  Loss 


500 

-0.02 

dB 

1000 

-0.04 

dB 

2000 

f'. 

o 

. 

o 

1 

dB 

4000 

-0.16 

dB 

A  loss  of  less  than  -0.2  dB  is  judged  to  be  quite  acceptable.  Thus, 
N'  =  4000  is  selected,  dividing  the  total  array  in  20  subarrays. 
(Actually,  it  is  the  random  phase  components  arising  from  element 
displacements  that  are  likely  to  limit  the  size  of  N' .  However, 
since  the  true  magnitude  of  such  displacements  remains  to  be  eval¬ 
uated,  N*  =  4000  is  assumed  for  the  ensuing  development.) 
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Returning,  now,  to  the  pulsed  wave  concept,  assume  that  a  pulse 
is  transmitted  whose  width  on  the  array  spans  4000  elements.  (Such 
a  pulse  would  have  a  very  practical  width  of  0.4  ys . )  Unlike  the 
very  narrow  pulse  whose  return  from  each  element  is  time  separable 
at  the  receiver  (Figs.  6  and  7) ,  the  return  from  4000  simultaneously 
illuminated  elements  produces  time  overlapping  signals  with  respect 
to  the  specified  subarrays.  Thus,  discrimination  between  the  sub¬ 
arrays  on  the  basis  of  simple  time  gating  is  not  possible.  Theoreti¬ 
cally,  a  means  is  needed  to  illuminate  a  4000-element  segment  of  the 
array,  and  then  instantly  hop  the  illumination  to  the  next  4000-ele¬ 
ment  segment,  etc.  One  way  of  accomplishing  such  illumination  could 
be  to  use  a  transmitting  aperture  of  sufficient  size  to  cause  a  beamwidth 
so  narrow  that  only  4000  elements  are  subtended  at  a  time.  It  should 
be  apparent  that  this  approach  is  more  unworkable  than  the  narrow  pulse 
method  because  it  would  require  an  aperture  larger  than  the  array  itself. 

The  practical  solution  is  based  on  the  use  of  the  correlation 

properties  of  pseudo— noise  (PN)  codes  [5].  Rather  than  transmitting  a 

pulse  whose  width  covers  4000  elements,  a  continuous  wave  carrier 

is  biphase  modulated  by  a  PN  code  whose  chip  period  spans  the  4000 

elements.  For  the  previously  cited  parameters,  a  PN  chip  rate  of 

2.5  MHz  is  required.  Furthermore,  because  modulation  of  a  continous 

wave  is  involved,  P  =  P  Thus,  the  pulse  width  and  power  speci— 

pulse  cw 

fications  are  well  within  attainable  limits. 

Consider  a  signal  of  the  form: 

s  (t)  =  PN ( t)  cos[u  t]  ,  ((>0) 

T  o 

input  to  the  array  model  of  Fig.  2.  From  Eq.  (15),  the  received  sig¬ 
nal  will  be: 

N-l 

Sp(t)  =  PN(t  -  nx)  cos[V  +  ^n+]J  > 

n=0 

where  x  =  T  +  T 2,  and  the  PN  chip  period  is  N'x.  If  the  received 
signal  given  by  (61)  is  correlated  at  the  receiver  with  PN(t),  the 


(61) 
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result  ing  correlation  function  obtained  is  (the  average  being  taken 
only  on  the  PN  components)  : 

r  (nx;  t)  =<s_(t)  PN(t)> 
o  R 

N-l 

=  ^  <PN(t)  PN (t  -  nx ) >cos  [V  +  W 
n=0 
N-l 

-  L  rPM(nT)  COSl"oC  +  *n+l!  •  ( 

n=0 


The  term  r  (nx)  is  the  autocorrelation  of  the  PN  code,  and  is  shown 
PN 

in  Fig.  8.  In  effect,  this  function  weights  the  cosine  terms  of  Eq. 
(62)  so  that  Eq.  (62)  may  be  written  as: 


r 

o 


(nx; 


N# 

n=0 


1  - 


n 

N' 


cos  [go  t  +  <j>  ,  ,  ]  . 
o  n+1 


(63) 


It  is  readily  seen  that  (63)  involves  the  returns  from  only  the  first 
Ny  =  4000  elements. 

Next,  correlate  (61)  with  PN(t  -  N*x)  .  The  result  becomes: 


N-l 

r^nx;  t)  =  £  rpN(N-x)  <nT>  cos  [^t  +  4>n+11  .  (64> 

n=0 

The  term  rpN(N/T)  (nT)  is  drawn  in  Fi8*  9,  where  rpN  (nx)  is  also 
shown  for  comparison.  By  inspection  of  Fig.  9,  (64)  is  written  as: 

N' 

rpnx;  t)  -  £  W  cos[V  +  Vl1 

n=0 

2N' 

+  E  4  -  257)  cos[v +  Vi1  •  (65) 

n=N' 

* 

The  autocorrelation  has  been  normalized  to  its  maximum  value. 


-31- 


Fig.  9  — 
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It  can  be  seen  that  (65)  involves  the  returns  from  the  first  2N#  = 

8000  elements.  Note,  however,  that  the  weightings  over  the  first 
N'  returns  are  just  the  opposite  from  the  weightings  in  Eq .  (63). 

In  a  similar  fashion,  additional  correlations  are  made  at  multi¬ 
ples  of  N't.  Each  additional  correlation  will  produce  a  weighting 
triangle  like  that  of  Fig.  9,  but  displaced  in  NT  by  the  multiple 
of  N't.  Thus,  each  correlator  (except  for  the  two  end  correlations) 
produces  a  weighted  sum  of  the  dephased  returns  from  2N'  elements. 
Because  the  adjacent  correlation  functions  overlap,  the  correlator 
outputs  are  not  independent  since  N'  elements  are  in  common.  Never¬ 
theless,  sufficient  discrimination  is  obtained  so  that  the  output  of 
any  given  correlator  can  be  represented  in  the  form: 

s  (t)  =  A  cos  [to  t  -  ip  ]  ,  (66) 

rC  K  OK. 

N 

k=D1  2  •••  —— 

K  UJ-L>Z>  >  N/ 

The  amplitude,  A^,  and  phase,  are  both  deterministic  functions 

of  the  quadratic  phase  terms,  <j)  ,  and  the  weighting  triangle.  Thus, 
may  be  corrected  relative  to  a  specified  reference.  Figure  10 
illustrates  a  functional  PN  correlation  processing  receiver  that  per¬ 
forms  all  of  the  above  indicated  operations.  The  bandpass  filters 
(BPF)  effectively  perform  the  correlation  averaging  or  weighting. 

Each  channel  is  phase  corrected  to  the  value  . 

A  practical  PACSAT  processing  receiver  must  perform  several 
functions.  In  addition  to  making  the  needed  near-field  phase  cor¬ 
rections,  it  must  also  phase-lock  to  the  carrier  and  demodulate  the 
message  data  placed  on  the  carrier.  (Data  are  assumed  to  be  biphase 
modulated  onto  the  cw  carrier  in  a  manner  identical  to  that  of  the 
PN  modulation.)  Figure  11  shows  a  practical  method  for  accomplishing 
all  three  functions.  Only  a  single  channel  of  the  (N  +  1)/N/  channels  is 
depicted  in  detail,  with  the  interfaces  to  all  other  channels  as  indi¬ 
cated. 


Practical  PACSAT  processing  receiver 
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PN  correlation  and  averaging  takes  place  in  a  combination  of 
bandpass  and  lowpass  filtering.  Following  the  PN  correlator,  the 
signal  is  demodulated  in  the  usual  inphase  ( 1^)  and  quadrature  phase 
(Q  )  manner.  Subsequent  to  demodulation  and  lowpass  filtering  (LPF) , 
phase  correction  is  implemented  by  means  of  arithmetic  operations  on 
the  I  and  Q  components,  as  shown  in  Fig.  11.  The  phase  correction 

1C  K. 

involves  the  functions: 


1  =  1  cos  ik  +  Q  sin  ip  (67) 

k  k  k  k 

and 

Q  =  I  sin  ip  -  Q  cos  ip  .  (68) 

k.  k  Ic  k 

Note  that  no  physical  phase  shifters  are  used;  the  processing  will 
most  likely  be  carried  out  using  digital  circuits. 

Carrier  coherence  is  obtained  by  means  of  a  Costas-loop  config¬ 
uration.  The  product  IQ  forms  the  carrier  tracking  error  function, 
which  is  summed  with  similar  IQ  terms  from  the  other  channels.  The 
composite  tracking  error  is  then  filtered  and  applied  to  the  input 
of  a  voltage  controlled  oscillator  (VCO) .  Output  of  the  VCO  is  phase- 
split  into  quadrature  components  and  applied  to  the  1^  and  Q^  demodu¬ 
lators,  thus  closing  the  tracking  loop. 

Finally,  the  I  terms  from  each  channel  are  summed  to  form  the 
demodulated  data  component,  which  is  subsequently  symbol  synchronized 


and  detected. 
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V.  CONCLUSIONS  AND  FURTHER  WORK 


This  Note  has  examined  the  problems  and  some  possible  solutions 
to  the  signal  acquisition  and  receiver  processing  of  PACSAT  relayed 
messages.  Although  much  groundwork  has  been  laid,  more  effort  is 
needed  to  fully  analyze  and  develop  the  concepts.  In  particular, 
more  work  should  be  accomplished  in  the  following  areas. 

1.  Using  the  PN  method  of  subarray  signal  discrimination,  a 
complete  analysis  of  the  theoretical  gains  obtainable  as 
a  function  of  subarray  size  and  PN  correlation  dependence 
is  needed. 

2.  Although  the  likelihood  of  random  phase  errors  due  to  arra 
element  displacement  has  been  mentioned,  no  specific  solu¬ 
tion  for  possible  correction  has  been  proposed.  Thus,  an 
investigation  of  the  means  and  degree  to  which  the  random 
components  may  be  estimated  and  corrected  is  required. 

3.  A  complete  functional  development  of  a  total  processing 
receiver  should  be  undertaken.  Aspects  to  be  included  are 

a.  Means  for  PN  code  synchronization. 

b.  Frequency  acquisition/detection  details. 

c.  Automatic  gain  control. 

d.  Digital  processing  algorithms. 

Implementation  risk  and  cost  estimates. 


e. 
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Appendix 

AN  EFFICIENT  CIRCUIT  FOR  FINDING  A  CARRIER  IN  A  WIDE  BANDWIDTH^ 

I.  S.  Reed 


INTRODUCTION 

Some  time  ago  the  author  of  this  appendix  introduced  a  cir¬ 
cuit  for  detecting  a  sine  wave  of  unknown  frequency  in  a  wide 
bandwidth  [l,2,3].  The  circuit  was  called  an  autocovariance  fil¬ 
ter.  The  detection  criterion  associated  with  this  circuit  has 
been  termed  semicoherent  detection. 

The  purpose  of  this  appendix  is  first  to  suggest  that  this  cir¬ 
cuit  or  some  variant  of  it  might  be  used  in  a  PACSAT-like  system  to 
efficiently  detect  and  locate  the  carrier  frequency  of  a  potential 
sender  in  order  to  establish  a  communication  link.  The  bandwidth  in 
which  such  a  carrier  can  be  is  quite  wide.  Also  the  input  signal-to- 
noise  ratio  of  the  carrier  is  only  moderate.  Hence  a  second  purpose 
is  to  compute  an  estimate  of  the  expected  output  signal-to-noise  ratio 
for  this  detection  criterion. 


SEMICOHERENT  DETECTION  AND  LOCATION  OF  A  CARRIER 

To  minimize  the  mathematical  description  of  the  semicoherent  de¬ 
tection  criterion,  assume  the  unknown  carrier  signal  is  complex  of  form 


s(t) 


,  i27rft  .  io)t 
Ae  =  Ae 


(A.  1) 


where  co  =  27rf  and  f,  the  unknown  carrier  frequency,  lies  in  the  band¬ 
width  interval. 


|  <  f  <  |  (A. 2) 

Physically,  signal  s(t)  is  obtained  first  by  heterodyning  a  real  sine 
wave  Asin(2tft)  down  in  frequency  so  that  the  band  of  interest  of  band¬ 
width  B  is  centered  about  zero  frequency,  as  given  by  Eq.  (A. 2).  The 

_ 

The  symbols  used  in  this  appendix  do  not  necessarily  conform 
with  those  used  in  the  text. 
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imaginary  part  of  Eq .  (A. 1)  is  obtained  as  the  quadrature  component  in 
the  heterodyning  process. 

Accompanying  any  signal  process  must  be  a  noise  process.  Origi¬ 
nally  this  noise  process  is  real  but  if  the  above  heterodyning  process 
is  applied  to  the  sum  of  the  signal  and  the  noise  processes,  one  ob¬ 
tains  the  complex  signal  plus  noise  process, 

w(t)  =  s(t)  +  z(t)  =  Ae”*^  +  z  (t)  ,  (A. 3) 

where  now  both  the  signal  Ae1*"^  and  the  noise  process  z(t)  lie  in  the 
bandwidth  interval  (Eq.  (A.2)). 

Assume  the  complex  noise  process 

z(t)  =  x ( t )  +  iy (t) 

is  the  classical  zero-mean  Gaussian  complex  process  where  y(t)  is  the 
Hilbert  transform  of  x(t).  This  is  the  usual  process  selected  for  an 
analysis  of  a  communication  channel.  Also  let 


R(t)  =  E  j^z  ( t )  z  (t  -  T)  J  , 

(where  E  denotes  expected  value  operator  and  the  bar  denotes  complex 
conjugate)  be  the  covariance  of  the  process  z(t)  with  itself,  some¬ 
times  called  the  autocovariance  of  z(t).  It  is  well  known  that  R(t) 
is  the  inverse  Fourier  transform  of  the  spectral  density  G(f)  of  the 
process  z(t),  i.e.. 


-CO 

R(T)  =  f  e27TifT  G  (f )  df 

J  —  oo 


With  these  preliminaries,  the  semicoherent  detection  criterion  for 
detecting  the  signal  Ae10Jt  in  noise  z(t)  is 


w(t)  w(t  -  t  )  dt 
o 


>  const 


(A. 4) 


-41- 


where  w(t)  is  given  by  Eq.  (A. 3),  T  is  called  the  integration  time, 

and  x  is  a  constant  delay  time  to  be  chosen.  If  t  =0,  then  Eq. 
o  o 

(A. 4)  becomes  the  classical  energy  detection  criterion. 


2 

r  (t)  dt  >_  const  (A. 5) 


where  r(t)  is  the  envelope  process  of  w(t).  The  disadvantage  of  Eq. 
(A. 5)  over  Eq.  (A. 4)  is  that  all  knowledge  of  the  phase  of  the  signal 
has  been  lost  in  the  statistic 


I 


y  =  I  w(t)  w(t  -  T  )  dt 


(A. 6) 


A  better  choice  for  T  is  to  choose  T  so  that  the  autocovariance 

o  o 


of  z(t)  is  zero,  i.e., 


R(Tq)  =  0, 


(A. 7) 


or  very  nearly  so.  Taking  the  expected  value  of  complex  statistic  y 
in  Eq.  (A. 6)  yields 


E(y)  =  f  E^w(t)  w(t  -  TQ)Jdt 

f  E^Ae1^  +  z(t)J  £ 

Jr\ 


,  ,  -ia)(t-T  ) 

Aelt0t  +  z(t)  Ae  °  +  z( 


1  -  r0>l 


■  f «[ 


n  IWT  -lto(t-T_)  .  .—7 - - \  loot  .  \~~7~z - 

2  o  ,  .  .  °  +  Az(t  -  T  )e  +  z(t)z(t  -  x  ) 

A  e  +  Az (t)e  o  o 


id  t 
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Using  Eq.  (A. 7)  and  the  fact  that  z(t)  is  zero  mean,  i.e., 

E[Z(t)]  =  E[ZU  “  Tq)]  =  °’ 


one  obtains 


iuvr  2  27TifT0 

E (y)  =  AZTe  °  =  A  Te 


(A. 8) 


as  the  mean  or  expected  value  of  statistic  y.  Observe  in  Eq.  (A. 8) 
that  the  frequency  of  the  carrier  can  be  obtained  from  the  complex 
number  E(y)  by  measuring  the  angle  2TrfTQ.  Similarly,  statistic  y  can 
be  used  to  find  an  estimate  of  f.  One  such  estimator  for  f  is 


(A. 9) 


where  Re(y)  and  Im(y)  are  the  real  and  imaginary  parts  of  y. 

The  two  hypotheses  of  the  detection  criterion  in  Eq.  (A. 4)  are 


H 


o: 


w(t) 


H^:  w(t) 


z ( t ) ,  the  noise-only  hypothesis 

4-  z(t),  the  signal-plus-noise 


hypothesis . 


The  semicoherent  detection  process  in  Eq.  (A. 4)  uses  a  nonlinear 
statistic.  Thus  for  the  figure  of  merit  of  such  a  nonlinear  statistic 
one  can  use  an  output  signal-to-noise  ratio  of  the  form 


(s/N)o 


lyy)  -  yy)i2  _  iyy>  -  y^i2 

E0|y-E0(y)|2  Eq|  y  | 2  -  |  EQ(y)  I  2 


(A. 10) 


where  EQ (y)  and  (y)  denote  the  conditional  expectations 
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Eq (y)  =  E[y  j  Hq]  and 
E-^y)  =  E[y  |  Hx ]  , 

with  respect  to  hypotheses  and  respectively. 

To  compute  Eq.  (A. 10)  for  statistic  y  note  first  that 

EQ(y)  =  0  (A. 11) 


and  by  Eq.  (A. 6)  that 


?  ia)Tn 

Ex(y)  =  A  Te 


(A. 12) 


Hence 


E0|y  -  E0(y)l2  =  E0lyl2  (a. 13) 

To  find  EQ | y |  in  Eq.  (A. 13)  in  a  simple  manner  it  is  desirable 
to  use  some  approximations.  Let  the  integral  for  y  in  Eq.  (A. 6)  be 
approximated  by  a  Riemann  sum  as  follows 


y 


Z ( t  —  T  )  dl 
O 


n 

z(kA)  z(k  -  p)A 

k=l 


(A. 14) 


where  A  is  the  reciprocal  of  the  bandwidth  B,  i.e.. 


A 


I 

B 


(A. 15) 


T  =  nA 
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and  p  is  the  least  integer  greater  than 


T 

O 


A 


T  B 
o 


(A. 16) 


P  =  [tqb] 

the  so-called  principal  part  of  number  T0B.  Let  the  spectrum  of  z(t) 
be  zero  outside  the  interval  — B  <_  f  <_  B.  That  is,  assume  the  spectrum 
is  band-limited.  Then 


R(A)  =  R(l/B)  =  0 


Letting  z(kA)  =  z^»  then  Eq.  (A. 17)  yields 

fe]  -  2°2Aj 

where 

a 2  =  E(Re  zjp  =  E  (Im  z^) 
and  6  is  the  Kroneker  delta  function  with 

ij 


(A. 17) 


(A. 18) 


Putting  the  sum  approximation  (Eq.  (A. 14))  in  Eq .  (A. 13)  produces 
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n  n 


e0 I y 1 2  -  A2 


'jL-J  E  (ZkZk-p  ZjZj-p) 


(A. 19) 


k=l  j=l 


Since  the  are  jointly  Gaussian,  one  may  use  the  identity  [A. 4], 


E(z^z^z^z^)  =  E  (z]_z2)  E  (Z3Z^)  +  e(ziz4^  e(z3Z2^’ 


on  Eq.  (A. 19),  thereby  obtaining 


E0 1 y 1 2  -  A2 


+  E 


£Z|'(v.-,) 

j=l  k=l  1 

(vj)^  ZJ-J| 


(A. 20) 


This  latter  expression  is  readily  evaluated  with  Eq.  (A. 18)  to  yield 
finally 


E0 1 y i 2  -  A2 


k=l X  '  N  / 


(A. 21) 


If  one  substitutes  Eqs.  (A. 11,  A. 12,  and  A. 21), 


2  iwT  2  /  2  \2 

(S/N)o  -  |At^  2°2  -(72)  !  ■  (S/N)i  BT 

AT (2a  )  \2a  / 


(A. 22) 


as  the  output  signal-to-noise  ratio  in  terms  of  the  usual  input  signal- 


(S/N)  =  ~ 

2a 


to-noise 
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This  result  is  the  same  figure  of  merit  measure  that  is  usually  obtained 
for  the  classical  square  law  or  power  detector  criteria  [5]. 

SOME  ADDITIONAL  REMARKS 

At  the  suggestion  of  the  author  of  this  appendix,  the  semicoherent 
criterion  or  autocovariance  was  tested  and  applied  to  long  pulse  de¬ 
tection  by  John  L.  Allen  [5]  at  the  MIT  Lincoln  Laboratory  in  1959. 

John  Allen's  circuit  was  attached  to  the  receiver  of  the  400  megahertz 
experimental  satellite  detection  radar  at  Millstone  Hill,  Massachusetts. 

The  Millstone  radar  had  a  long  cw  pulse  about  one  millisecond 
in  duration.  The  doppler  filter  bank  to  cover  the  doppler  band  of  in¬ 
terest  required  about  512  analog  filters,  overlapping  at  the  three 
decibel  level,  with  each  filter  matched  to  a  one  millisecond  pulse. 

The  semicoherent  detection  circuit  was  put  in  parallel  with  this  filter 
bank  for  comparative  evaluation. 

For  the  same  measured  false-alarm  rate,  the  autocovariance  detec¬ 
tion  and  doppler  estimation  circuit  appeared  to  perform  almost  as  well 
as  the  doppler  filter  bank  system  in  terms  of  probability  of  satellite 
detection.  Every  time  the  filter  bank  detected  a  satellite,  John 
Allen's  circuit  also  made  a  detection. 

Other  more  pressing  matters  and  the  vicissitudes  of  the  time  made 
it  difficult  to  accomplish  carefully  controlled  experiments  with  semi¬ 
coherent  detection  on  the  Millstone  Hill  radar.  For  one  thing,  the 
increased  number  of  satellites  made  this  radar  an  operational  radar, 
rather  than  an  experimental  tool  of  the  laboratory.  Finally,  because 
of  the  analog  nature  and  complexity  of  the  filter  bank,  its  performance 
would  vary  almost  on  a  daily  basis,  making  any  comparison  with  another 
detection  criterion  uncertain  at  best.  In  spite  of  these  deficiencies 
in  the  experiment,  the  informal  results  seem  to  indicate  that  semi¬ 
coherent  detection  is  qualitatively  comparable  in  terms  of  probability 
of  detection  with  a  fixed  false-alarm  rate  to  square  law  detection. 

The  semicoherent  detector  can  be  implemented  at  the  complex  video 
level,  or  at  the  output  of  the  I.F.  section  of  the  radar  receiver. 

John  Allen's  original  circuit  was  realized  by  the  latter  approach  using 
mixers  and  analog  delay  lines  and  integrators.  A  simplified  block 
diagram  of  this  circuit  is  shown  in  Fig.  A.l. 


Fig.  A-1  — Simplified  block  diagram  of  a  semicoherent  detector  operating  on  an  I.F.  waveform 
( the  phases  of  noise-free  signals  are  indicated  at  different  points  in  the  diagram ) 
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